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THE HEIGHT OF MULTIPLE EDGE PLANE TREES 


CLEMENS HEUBERGER, HELMUT PRODINGER, 
AND STEPHAN WAGNER 


Abstract. Multi-edge trees as introduced in a recent paper of 
Dziemianczuk are plane trees where multiple edges are allowed. 
We first show that d-ary multi-edge trees where the out-degrees 
are bounded by d are in bijection with classical d-ary trees. This 
allows us to analyse parameters such as the height. 

The main part of this paper is concerned with multi-edge trees 
counted by their number of edges. The distribution of the number 
of vertices as well as the height are analysed asymptotically. 


1. Introduction 

Dziemianczuk [4] has introduced a tree model based on plane (=pla- 
nar) trees m p. 31], which are enumerated by Catalan numbers. In¬ 
stead of connecting two vertices by one edge, in his multi-edge model, 
two vertices can be connected by several edges. If one counts trees 
by vertices, one must somehow restrict the number of edges in order 
to avoid an infinity of objects with the same number of vertices. In 
[4] , the chosen restriction is that each vertex has out-degree at most d, 
i.e., there are at most d edges going out from any vertex. However, if 
one counts trees with a given number of edges, the restriction with the 
parameter d is no longer necessary. This is in contrast to the case of 
classical plane trees where the number of edges equals the number of 
vertices minus one. 

In [4], several parameters of multi-edge trees were analysed, but some 
questions about the (average) height (i.e., the maximum distance from 
the root) of such multi-edge trees were left open. The present paper 
aims to close this gap. 
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In Section [2j a Injection is constructed which links d-ary multiple 
edge trees with standard d-ary trees. Since the bijection is height¬ 
preserving, and the height of d-ary trees is well understood, we can 
resort to results by Flajolet and Odlyzko [7] as well as by Flajolet, 
Gao, Odlyzko and Richmond £5] and provide in this way a full analysis 
of the height of d-ary multi-edge trees, cf. Theorem |2.3 


In Section [3j we count trees by the number of edges and drop the 
parameter d. The analysis of the height of plane trees appears in a 
classic paper by de Bruijn, Knuth and Rice [1] (see also [H]), with 
an average height of asymptotically y/wn. Now, we can follow this 
approach to some extent, but combine it with a technique presented 
in [9j. The expected height is asymptotically equal to y/nn, with 
a more precise result in Theorem |3.6[ The constant is smaller, which 


is also intuitive, since the multiple edges contribute to the size of the 
objects, but not to the height. We also give an exact counting formula 
in terms of weighted trinomial coefficients (Theorem 3.5) and a local 
limit theorem (Theorem 3.8). 

The distribution of the number of vertices in plane multi-edge trees 
with n edges is analysed in Theorem 3.11 The number of trees with 


given number of vertices and edges is given in Theorem 3.10 


2. A BIJECTION BETWEEN d-ARY MULTI-EDGE TREES AND 
ORDINARY d-ARY TREES 

As explained in the introduction, Dzicmianczuk [3J studies d-ary 
multi-edge trees, where a vertex can have at most d edges going out 
from it. We present a simple bijection to ordinary (pruned) d-ary trees, 
where every vertex has d possible positions for an edge to be attached 
(e.g., left, middle, right in the case d = 3). See [10j, Example 1.14] for 
a discussion of pruned d-ary trees. This bijection preserves (amongst 
other parameters, such as the number of leaves) the height, allowing us 
to reduce the problem of enumerating d-ary multi-edge trees by height 
to the analogous question for d-ary trees, which has been settled in jS], 

Our bijection can be described as follows: suppose that a vertex v 
of a d-ary multi-edge tree has r children, which are connected to v by 
ki, k 2 , ..., k r edges respectively. The corresponding vertex v' in the 
d-ary tree also has r children (corresponding to the children of v in the 
natural way), which are attached to v' by edges in the Aq-th, (£q + k 2 )- 
th, (ki + k- 2 + k 3 )-th, ..., (ki + k 2 + • ■ ■ + k r )-th position. Since we are 
assuming that ki + k 2 + • • • + k r is always < d, this is possible, and 
clearly this process is bijective for each vertex, so it also describes a 
bijection between trees. Figures [l] and [2] illustrate an example in the 
case d = 5. 

From this bijection, we immediately obtain the following corollaries: 
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Figure 1. A 5-ary multi-edge tree. 



Figure 2. The associated 5-ary tree, where each vertex 
can have children in five different positions (far left, left, 
middle, right, far right). 

Corollary 2.1. The number of d-ary multi-edge trees with n vertices 
equals the number of d-ary trees with n vertices, which is the Fuss- 
Catalan number - ( nd f]. 

n \n—lJ 

This is for instance shown in [10) Example 1.14], 

Corollary 2.2. The number of d-ary multi-edge trees of height h with 
n vertices equals the number of d-ary trees of height h with n vertices. 

It is well known that d-ary trees belong to the general class of simply 
generated families of trees, and the height of such families was studied 
in great detail in a paper by Flajolet, Gao, Odlyzko and Richmond {5] . 
They obtain the following local limit theorem (only stated for d-ary 
trees here, i.e. setting <f(y) = (1 + y) d and r — 1/(d— 1) in the formulae 
given there), which refines earlier results of Flajolet and Odlyzko [7] 
on the average height: 

Theorem 2.3 ([5j Theorem 1.2]). Let Njf\n) be the number of d- 
ary trees (d-ary multi-edge trees) with n vertices whose height is h and 
N^ d \n) the total number of d-ary trees (d-ary multi-edge trees) with n 
vertices. For any S > 0, we have the asymptotic formula 



2c/3 4 \Ztt/ 7i'y^ j (7mr) 2 (2 (tt mf3) 2 — 3) e ( 


(7rm/3) 2 


NW{n) 


m> 1 


2c/r 1 n~ 1/2 ^ m 2 (2(m/P) 2 - 3)e " (m/ ^ 


m> 1 
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where c = a/ 2 (d — 1 )/d, uniformly for 

d'/logn) -1 / 2 < (3 = 2 y/n/(ch) < d(logn) 1 / 2 . 

Corollary 2.4 ([71 Theorem S], [5], Corollary 1.2]). The average height 
of d-ary trees with n vertices (and thus also the average height of d-ary 
multi-edge trees with n vertices) is asymptotically equal to 

\j2ndn/ (d — 1). 

Similar results for the average height were obtained by Kemp [131 
M\ (see also HZ!) for slightly different models of random plane trees, 
namely for trees with given root degree or number of leaves. 

As it was mentioned earlier, other statistical results carry over from 
d-ary trees to d-ary multi-edge trees as well: 

Corollary 2.5. The number of d-ary multi-edge trees with n vertices 
and k leaves equals the number of d-ary trees with n vertices and k 
leaves. 

More generally, the following holds: 

Corollary 2.6. For every r € {0,1,..., d}, the number of d-ary multi¬ 
edge trees with n vertices, k of which have exactly r children, equals the 
number of d-ary trees withn vertices, k of which have exactly r children. 
Thus the average number of vertices with exactly r children is the same 
for d-ary multi-edge trees and d-ary trees with n vertices. 

It is not difficult to show that the average proportion of vertices with 
exactly r children is asymptotically equal to (d— l) d ~ r d~ d as n —» oo 
(cf. the paragraph following [3], Theorem 3.13] with k — r, (fk — (,), 
$(y) = (1 + y) d , r = l/(d — 1)), which tends to l/(r!e) as d — y oo. 
This generalises the observation made in [4] in the case r = 0 that the 
asymptotic average proportion of leaves tends to 1/e as d —> oo. 


3. Trees with Given Number of Edges 


In this section, we consider plane rooted multi-edge trees with a given 
number n of edges (which we call the size of a tree). The resulting 
counting sequence A n is sequence A002212 in [15], see also [|T6]. It 
starts with 1, 1, 3, 10, 36, 137, 543, 2219, 9285, 39587. 

Asymptotically, the number A n of plane rooted multi-edge trees with 
n edges is 


( 1 ) 


A r 


^L(i + oO-)\ 

2^/irn 3 \ \nJJ 


This will follow without further effort at the end of the proof of Theo¬ 
rem 


3.6 We now analyse the height of multi-edge trees. 
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3.1. Generating Functions. In the following lemma, we introduce 
the fundamental transformation which will be used throughout this 
section. The principal branch of the square root function is chosen as 
usual, i.e., as a holomorphic function on C \ M< 0 such that \/l = 1. 


Lemma 3.1. Let Z — C \ [1/5,1] and U = {u 6 
(—3 + \/5)/2}. Let 

1 — 32 — a/1 — 5z\/l — z 


V[Z = 


C(«) = 


2z 


u 


u 2 + 3u + 1 


for z G C, 
for u G C \ 


u\ < l;u / 


3± 

J 


T/ien v: Z ^ U and (: U —>■ Z are bijective holomorphic functions 
which are inverses of each other. 


Proof. We first note that is well-defined and holomorphic on U with 
C'( u ) V 0 for all u G U. If |w| = 1, then 

, 1 1 
^ U) u+ l + 3 3 + 25 lu 

U 

Thus the image of the unit circle under ( is the interval [1/5,1]. 

For every z G C \ {0}, z = f{u) is equivalent to 

( 2 ) u 2 + u( 3—^)+1 = 0 

which has two not necessarily distinct solutions u±, u-2 G C with U\U2 = 
1. W.l.o.g., |wi| < \uf\. Thus either u\ G U and |i*21 > 1 or |«i| = 
\v/ 2 1 = 1. In the latter case, we have z G [1/5,1]. For z = 0, z = 
f{u) is equivalent to u — 0. This implies that (: U -A Z is bijective. 
Furthermore, (: U —^ Z has a holomorphic inverse </ _1 dehned on the 
simply connected region C \ [1/5, oo). 

Solving (j2]) explicitly yields 

1 - 3z ± Vl-6z + 9z 2 - 4z 2 1 - 3z ± y/1 - 5 Zy/l^Z 

u - 2z “ 2 ~z ' 

In a neighbourhood of zero, we must have C _1 (^) = v ( z )> because 

1 - 3^ + \/l - 5z^/l^z 
2 z 

has a pole at z = 0. 

It is easily seen that V1 — 5z^/l — z is a holomorphic function on Z. 
By the identity theorem, (V 1 = v holds in C \ [1/5, oo). By continuity 
of v in Z, v is also the inverse of ( in (1, oo). □ 

For h > 0, consider the class 7/ of plane rooted multi-edge trees of 
height at most h. Denote the ordinary generating function associated 
to T h by T h {z). 
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Lemma 3.2. The generating function Th(z) is given by 


( 3 ) 

where 

( 4 ) 


T h (z) = (1 -z) 


a 


a 


= (« + 1)t- 


u 


h +1 


h+2 


P h +2 


U 


h+2 


a = 


1 — z + \/l — 5z^l — z 


u + 1 


P = 

for z = C (u) G Z. 


2 u 2 + 3u + 1 ’ 

\/l — 5z^l — z u{u + 1) 

2 u 2 + 3u + 1 


Proof. The class % consists of an isolated vertex. For h > 0,Th consists 
of a root and a sequence of branches of height at most h — 1 such that 
each branch is attached by a positive number of edges to the root. If 
£ = {e} is the class of one edge, we can write Th symbolically as 

(5) T h = ox (£+%-!)*. 

The symbolic equation ([5]) translates to 

1 1 — z 

h * 1 - j^T h _i(z) 1- z - zT h _i(z )' 

This may be seen as a continued fraction. To obtain an explicit 
expression for Th(z), we use the ansatz Th(z) = Ph{z)/qh(z ) with 
Po(z) = qo{z ) = 1 and 

p h (z) = (1 - z)q h - l (z), 

q h (z) = (1 - z)q h _ l (z) - zp h -i(z). 

Eliminating Ph{z) yields the second order recurrence 
q h (z) = (1 - z)q h _ 1 (z) - z( 1 - z)q h _ 2 (z). 

The characteristic equation is 

Q 2 - (1 - z)Q + z(l - z) = 0. 

This quadratic equation has the roots a and (3 defined in Q. This 
yields the explicit expressions 


Qh(z) = 

which result in 

( 6 ) 


a 


h+2 




h+2 


(1 - z)(a- py 


Th(z) = (1 - z) 


Ph(z) = - 

h -\-1 _ 

U l o 

a 


h+1 


-P 


h -\-1 


a — P 


a 


h+2 


P h +r 


Under the substitution z — ((u), we have 


z = 


(u + lf 


u 2 + 3u + 1 ’ 
Inserting this in Q yields ([3]). 


P = 


u{u + 1) 
u 2 + 3u + 1 ’ 


u+1 


a = 


u 2 + 3u + 1 


□ 
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Let T be the generating function of all plane, rooted multi-edge trees. 
Lemma 3.3. For z = ((u) G Z, 

T(z) = -= u+l 

z 

and 


(7) 


C T-T h )(z ) 


1 -u 2 u h+2 
u 1 — u h+2 


Proof. It is clear that T is the limit of T h for h -A oo. As |tt| < 1, we 
have T(z) = u + 1. The expression for T — Th follows. □ 


Note that T could also have been determined by removing the restric¬ 
tion on h in the symbolic equation and solving the resulting quadratic 
equation for T. 

Lemma 3.4. The functions T(z), T h (z ) and J2 h>0 (T — T h )(z) are 
analytic for z G Z. 

Proof. By the explicit formula for fd. it is clear that T(z) is an analytic 
function on Z. 

For u = v(z) and z G Z, the function 

1 — 

W = (“ + 1 ) r ^^ 


is clearly analytic. 

The sum ^2 h>0 (T — T h )(z) can be written as 




1 - w 2 


u 


E 

h> 0 


U h+2 

1 — u h+2 


We can bound the sum by 


E 


u h + 2 
1 — u h+2 



E l \h+2 

M 

h> 0 


(i - M 2 )(i - H)' 


By the Weierstrass M-test, 


Et-f.) w 

h> 0 

converges uniformly on compact subsets of U and is therefore analytic 
in U. 

The results for z G Z follow by the fact that v(z) is analytic. □ 
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3.2. Explicit Formula for the Number of Trees of Given Height. 

At this stage, we can compute the number of rooted plane multi-edge 
trees of size n and height > h explicitly. Taking the difference for h 
and h — 1 results in a formula for the number of trees of height h. 

Theorem 3.5. Let h > 0. The number of rooted plane multi-edge trees 
of size n and height > h is 


( 8 ) 


Ei 

fc >0 


where 


n - 1; 1, 3,1 
n — {h + 1) — (h + 2 )k 


- 2 


n - 1; 1, 3,1 

n — (h + 1) — (h + 2 )k — 2 
n - 1; 1,3,1 

n — {h + 1) — (h + 2 )k — 4 


71 ] 1, 3,1 
k 


2\n 


= [n fe ](l + 3n + n 2 ) 


denotes a weighted trinomial coefficient. 

Proof. By the definition of the generating functions, we have to com¬ 
pute [ z n ){T — T h )(z). By Cauchy’s formula, we have 


( 9 ) 


M(T-7W(*) = - 


(T - T h )(z) 


7 n+l 


dz. 


I 7' 1 

, \z\ small ^ 

For sufficiently small u|, the index of 0 with respect to f{u) is 1. 
Therefore, using the substitution z = ((u) and using Cauchy’s formula 
again, we can rewrite ([9]) as 


[z n }(T 


T h )(z) 

1 

2ni 


(T-T h )( CM) 


|u| small 


U 


n+1 


(u 2 + 3 u + l) n X (1 — u 2 ) du 


= [u n ] (T - T h )(((u))(u 2 + 3 u + l) n -\l - u 2 ) 

= K] ( 1 r_"y i ( “ 2 + 3 “+ir‘(i - “ 2 ). 

Expanding the denominator into a geometric series yields 
[z n ](T - T h ){z ) = [u n ] - u 2 ) 2 u h+lHh+2)k (u 2 + 3 u + l) n_1 


J2[u n ^ h+1+(h+2)k) }{ 1 - 2 u 2 + u 4 ){u 2 + 3 u + l) n_1 

fc >0 

Y. ( [ u n -( h + 1 +( h + 2 ) k l] ( u 2 + 3 u + l) n_1 


fc >0 


- 2[u n - (h+1+{h+2)k) - 2 ](u 2 + 3 u + l)”- 1 
+ [ u n -^ h+l +l h+2 ) k )- A ][u 2 + 3 u + l) n_1 ). 
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By the definition of 


n;l,3,l 


), this is exactly 


□ 


Remark 1. It would be possible to determine the asymptotic behaviour 
of the trinomial coefficients by means of the saddle point method (cf. 
dn Section 4.3.3]j and to obtain asymptotics for the average height 
(Theorem 3.6) and the local limit theorem (Theorem 3.8) from that, 
but the calculations would be somewhat more involved. 

3.3. Expected Height. We now compute the expected height of a 
random rooted plane multi-edge tree of size n. 

Theorem 3.6. Let H n be the height of a random rooted plane multi¬ 
edge tree of size n. Then 


e(j/„) = T^-? + o(T). 


( 10 ) 


Before proving Theorem |3.6[ we prove a lemma on the harmonic sum 
occurring in its proof. 


Lemma 3.7. We have 


(id £ 


u 


log(l 


u 


+ 


7 


h>l 


U“ 


U 


U 


+ log(i_^) _ l _ 2 + 0((1 _ u)log(1 _ tl)) 

as u —* 1 with |arg(l — u)\ < tt/ 3, where 7 is the Euler-Mascheroni 
constant. 


Proof. Using the substitution u = e 1 yields 


E li v ■'v 6 

1 — u h 1 — 


~—ht 


>—ht 


h> 1 


h> 1 


££ 

h> 1 k> 1 


g —kht _ 


d(m)e 


—mt 


m> 1 


where d(m) is the number of positive divisors of m. 
By [61 Example 11], we have 




—mt 


m> 1 


= i(— logt + 7) + 1 + 0 ( 1 ) 


for real t —> 0 + . However, the same argument can also be used for 
| arg 1 1 < 7 t /4 because the inverse Mcllin transform 


r»c+ioo 


-t 


= — r s T(s)ds 

‘ZTTZ .1 


remains valid for complex t with |argt| < 27r/5 by the identity theorem 
for analytic functions; cf. [9]. 

As 

(1 - u) 2 

t = - logu = - log(l - (1 - u)) = (!-«) + -— -h 0((1 - u ) 3 ), 


substituting back yields ( 11 ). 


□ 
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Proof of Theorem 3J3. We use the well-known identity 


E (H n ) = kHH n = k) = Y, P (#« = k ) = Yl P (^ > h ) 

k =0 k>h >0 h >0 


= - p (^» S h)) 

h> 0 


[z n ]T(z) 


We intend to compute [z 11 ] 'Yh h>Q {T — Th){z) via singularity analysis. 
The dominant singularity is at z — 1/5. To perform singularity analy¬ 
sis, we need the expansion of T — Th around z = 1/5, corresponding to 
u — 1 under the substitution z = ((u). 

By0 , we have 


]T(T-T h )K(u)) 

h> 0 


1 — u 1 


E 


u 


u —' 1 — u h 

h>2 


1 - U 2 U h 

-(! + «) + —]T — 


u 


h> 1 


U n 


IS 


By Lemma , this 

- T „)(<(«)) = —21og(l - u) - (2 - 2 7 ) 


( 12 ) 


/i>0 


+ 2 “ u ) + °(( 1 _ M ) 2 lo s ( 1 - “))• 


We have 


1 — u — — 5 z — -(1 — 5 z) + 0((1 — 5 z) 3 ^ 2 ), 

log(l-u) = ^log(l -5*) + ^log5- -5z + 0((l -5z)). 


Inserting this in (12) yields 


£(T-T fc )(z) = - log( 1 - 5~) - (2-2y +log5) 

h> 0 

+ ^V5VT^ + 0((1 - 5z) log(l - 5z)) 

for z -A | and |arg(| — z)\ < 37t/5, i.e. |arg(z — |)| > 2n/5. Note 
that the exact bounds for the arguments are somewhat arbitrary: the 
essential property of 27t/5 here is that it is less than 7t/2. Using the 
expansions of 1 —u and t in terms of y/1 — 5 z and of 1 — u, respectively, 
the angles are transformed accordingly, but we have to allow for a small 
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error. By singularity analysis [8], this yields 


3V5 5 

> (T-T h )(z) = — + 
n 


(13) 


h>0 


2 r(— l/2)n 3 / 2 


+ 0 5 


, log n 


n 


5 n 3 . 5 n+l/2 / logn 

- + 0 5 — 


n 


4a/- 


7in J 


n- 


The number of plane rooted multi-edge trees of size n is 
= [z n ]T(z) = [z n }(u + 1) = [z n }(2 - (1 - u)) 

= [z n } (2 - vWT^5z + |(1 — 5 z) + 0(( 1 - 5z) 3/2 )^j . 

Singularity analysis yields 


5 n n~ 3 / 2 f 5 

A n — —v5——- 7—7 + O 


r(— 1 / 2 ) 


n 5 / 2 ) 2 a/' 71 "'" 3 


5 n+ 1/2 + 0 n 

7rn° \ 


Combining this with (13) yields (10). 


□ 


3.4. Local Limit Theorem. In this section, we prove a local limit 
theorem for the height of a plane rooted multi-edge tree. As our gen¬ 
erating function is very explicit, we can give a result in a wider range 
than [5j. 


Theorem 3.8. Let 0 < e < |. Then, for 

Id' 

(14) 


4n7r 2 


< h < n 3/4-£ , 


5£ log n 

the probability of a plane rooted multi-edge tree to have height h is 

for 


G(a) = ^(2o 2 m 2 — 3)m 2 exp(— a 2 m 2 ) 


(15) 


m> 1 


cr 


m> 1 


5 ^ ( 2 0 2 ™ 2 - 3 ) “ 2 exp (- 0 “ 


The fact that the two expressions for G(ct) in (15) are equal is Pois¬ 
son’s sum formula (cf.|2j (3.12.1)]) for f(x) = (2 era?— 3)x 2 exp (—a 2 x 2 ). 

We first compute the integral which will appear by application of the 
saddle point method. 
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Lemma 3.9. Let 0<a<l, 0<6&e real numbers and c, d be complex 
numbers. Then 


dt 


(ct + d) 3 ex p(-f) 

(1 — ae lbt ) 2 

= + 1) ^cibm + d'jc 2 + Q cibm + d^j j 


m> 0 


x a m exp(-^ 2 m 2 ). 


Proof. We expand the denominator of the integrand as a binomial se¬ 
ries, dominated by (1 — a)~ 2 . Thus 



(ct + d) 3 ex p(-f) 
(1 — ae lbt ) 2 


= y(m + l)a m 

m> 0 


(cf + d) 3 exp 


f 2 

5 


+ ibmt ) 


dt. 


Substituting t = z + |ibm and shifting the path of integration back to 
the real line yields 


/ \ 

(ct + d) 3 exp — + ibmtj dt 

— ^6 2 m 2 j ^cibm + djc 2 + ^cibm + d)j j. 


= exp 


□ 


Proof. Instead of computing the number of trees of height exactly h , we 
compute the number A n h of trees of height exactly h — 1 because this 
leads to more convenient formulae and does not matter asymptotically. 
By 0, we get 


A nh := [z n }(T h _i - T h _ 2 ){z) = [z n ]((T - T h _ 2 ) - (T - T h _i)) 


u 


= z 


= Z" 


= Z" 


u 


u 


u 


h +1 


U 


- U" 

u h — u 


1 — u h+1 

h -\-1 


u (1 — u h )( 1 — u h+l ) 


(1 + w)(l - uf 


u 


u 


(1 — u h )( 1 — u h+l ) 


for z = ((u). 
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Using this transformation and Cauchy’s formula as in the proof of 

(1 + m)(1-m) 2 


Theorem 3J3 yields 
1 

2 m 
1 

2ni 


A-nh — 


M 


dz 

|*| small U (l-^)(l-^+ 1 )^+ 1 

(1 + m) 2 (1 - m) 3 (m 2 + 3 u + l) n_1 


W 


small (1 -«*)(! ~U h+1 )ti n - fc+2 


du. 


Now we apply the saddle point method to this integral. It turns out 
that the right choice for the contour of integration is given by the 
parametrisation u = re lip with r = exp(—“) and r < ip < 7r. This 
yields 


Anh ~ 2vr 


(1 + m) 2 (1-m) 3 (m 2 + 3m + 1) 

(1 — U h )(l — u h+1 ) U n -h+ l 

g(u) exp(nf(u)) 


\n— 1 


dip 


2n J_ n (1 — u h )( 1 — u h+l ) 


dip 


for 


f{u) = log(l + 3-U + M 2 ) + — lj log M, 

(1 + m) 2 (1 — m) 3 


g{v) = 


We set 


m(m 2 + 3m + 1) 

5 h 2 


a 2 = 


By the assumption (14), we have 
(16) 


4n 


7T 


a > 


e log n 


Note that r —* 1 for n —> oo. We also note that g[u) = 0(1) on the 
area of integration. If a 2 < 7r, 

|1 ~u h \ > 1 - r h = 1 -exp(-2 a 2 ) > 2a 2 exp(-2ct 2 ) > 2a 2 exp(-27r); 


thus 

(17) 

rhs = °(£) = 0(l0E ' re) ' = °(I0 = 0(logn) ' 

Otherwise, r h < exp(—27 t), i.e., ^0-^ and 3 h+1 are bounded. Thus 
can be used in any case. 

We first prune the tails. We set S n = n~ 1//2+e such that nS 2 = n 2e and 
n<5 4 = n~ 1+4e < n~ 1/,2+e and n§ n /h 2 = 0(mT 1//2+£ logn) for n —>• oo. In 
particular, we have 6 n = o(l). 
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For \(f \ > S n , we have 

|1 + 3u + u 2 \ < |1 + 3u| + r 2 = \/l + 6 r cos cp + 9r 2 + r 2 

< \/l + 6r cos <5 n + 9r 2 + r 2 

I £2 

< w 1 + 9r 2 + 6r — + r 2 = \/(l + 3r) 2 — 2r<5 2 + 


< 1 + 3r + r 2 — 


A 2 A2 

" < 1 + 3r + r 2 - 


1 + 3r 2 “ 10 

for sufficiently large n. We conclude that for \<p\ > 5 n , 

3R/(w) < log(l + 3r + r 2 - ^ - l) logr < /(r) - ^ 

for sufficiently large n. Thus, by 0, 


A n h — _ 


g(u)exp(nf(u)) 


dp 




2vr J_ Sn (1 - u h )( 1 - u h+1 ) 

+ C»(log 2 nexp(n/(r)) exp^— j. 
We now approximate the integrand in the central region. We have 

m = iog5+£(~+ip) + g(-|+^) 2 
+ o((^m)‘) 


5/r 2 




= 1 °g 5 - 73-ir + ° r + M 


= i+od + " H 


4n 2 5 

3 


h 


n 


M h 2 

Therefore, noting that n(h/n + <5 n ) 4 = 0(h 4 /n 3 + n<5 4 ) = o(l) yields 


A n h — 


2-5"exp(—f) f s - (f n -ivY 


57r 


!— s n ( i -« h ) 2 


exp 


h /i 4 log 


x (1 + O ( —I—, + 


n 


n n J 


n l ! 2 ~ e 
.2 


nip 




O 5" log n exp 


5/z 2 n<5 2 \ 

’ 4^" ~~ 100/ 
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We now use the substitution y/nip = t, leading to 

A nh 5n^ /^(i-^) 3 


2-5"exp(-^) J-s nV * (1 ~u h ) 2 


exp 


h h 4 log 


x 1 + 0 - + — + 


n 


n n J 


1 / 2 —£ 


n 


dt 


+ o(yn log 2 n exp ) • 


We set 


1hn 

Ehn 


r-oo (§L -i t \ 3 , 

-oo (1 -«0 2 d a) • 

5 h |t| W 


(1 ~r h ) 2 J_ 00 \2n 


- —I—7=) exp — — ) dt, 


and note that the contribution of \t\ > S n y/n is again negligible: we 
have 



(5h 
V2 n 



(1 — u h ) 2 


exp 



dt 


1 r°° /5 h 

- (1 -r h y Js Wn \2n + 




dt. 


Now we can use the estimate © for 1 — r h as before, and the inte¬ 
gral in the upper bound can in principle be computed explicitly. It is 
0((-y/n<5 n )“ 1 exp(—n<5 2 /5)), so the total contribution of the tails (i.e., 
the regions where \t\ > 5 n y/n; of course the estimate for negative t is 
analogous) is (^(n -1 / 2 ^ 1 log 2 nexp(— nS^/5)). It would be possible to 
give an even better bound, but this is enough for our purposes. 

We obtain 


(18) 


--'../.■I" \/if 
2-5"exp(- 5 £) 


= Ihn + 0 [E hn [ h - + K+ 108,1 


n rr 


1/2 —£ 


n 


+ O^Vn, log 2 n exp ) • 


By Lemma 3.9 with a = exp(—(5/r 2 )/(2n)), b = h/y/n , c = —i/yfn 


and d = (5 h) / (2 n) and by replacing m + 1 by m, we obtain 


( 19 ) 


25/ia/57t exp(^) ^/5 h 2 2 \ 2 

hn = - 4n/ 2^ ( IT' 171 ~ 3 ) m eX P 


4n 2 


25hy/5n exp(a 2 ) 


4n 2 


m> 1 


G{a). 


5 h 2 


4 n 


-m 




















16 


C. HEUBERGER, H. PRODINGER, AND S. WAGNER 


The integral Ehn can be bounded by 


( 20 ) 


Eh n ~ O 


h? j_1_ 

n 3 T n 3/2 


(1 


r h\2 


We first consider the case a 2 > tt. In this case, we have E hn = 
0(h 3 /n 3 ). All summands in the first expression in ( [15] ) are positive 
and its first summand is at least a 2 exp (—a 2 ), so that 

Ihn ~ ^(^ 3 ) = ^(Ehn)- 


Then (18) yields 

( 21 ) 

Anhh'Ky/n 25h\/57T exp(ct 2 ) 


h / 1 4 log n 


2-5-exp(-f r ) 


4n 2 


G(a) (l + o(— + —^ + 

\ \n n 6 


1 / 2 —£ 


n 


We now turn to the case a 2 < n. We now use the second expression 
for G(a) in (15). Again, all summands are positive and we bound G(a) 
by the first summand from below. This yields 

G(a) = nf^-expf-^- 


cr 


cr 


,3/2 


,3/2—e, 


and, by (19) and 

4n = ^(-^-exp(-elogn)J 

For an upper bound of Ehn, we use the estimate (1 — r h )~ l = 0(n/h 2 ), 
cf. ( Jl7[ ) . We get 


Ehn — o 


• SUo 


n 3 / 2 h 4 / 


77 , 


1/2 


/r 4 


■ E 3 / 2 e h 2 \ 

= 0, W^A ne ) =0( "' /to) - 


Thus (18) yields 


/ 00 n A nh bTi^n 25/iv / 5tt exp(ai 2 ) / / logn 

(22) 2.5-exp(-g) = - G <“)( 1+0 (^r 

Combining (21) and (22) with ([!]) yields the result. 


□ 


3.5. Number of Vertices. In this section, we consider the number of 
vertices of a random rooted plane multi-edge tree of size n. 

We first give an explicit formula. 


Theorem 3.10. The number of rooted plane multi-edge trees of size n 
with k vertices is 


( 23 ) 


1 (2k - 2\ f n - 1\ 
k\k-l )\k-2j' 
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Proof. We first provide a proof based on the generating function, which 
will also be needed later. Let T(y,z) be the bivariate generating func¬ 
tion for rooted plane multi-edge trees, where y marks the number of 
vertices and z the number of edges. Rooted plane multi-edge trees T 
can be represented symbolically as 

(24) T={y}x(S + T)*. 

This symbolic equation translates to 

y y( i - z) 


(25) 


T{y,z) = 


l-I ^ T (y,z) 1 -z-zT(y,z)' 

For a fixed z, we compute the coefficient [y k ]T(y, z) using the La¬ 
grange inversion formula. By (25), we have 

rp, si ~z-zT(y,z) 
y = T{y, z) -—-■ 

-L Z 

Now the Lagrange inversion formula gives us 

1 —!. i, / 1 — z 




zT 


= k [Tk ~ 1] 


1 - -+-T 

1—Z 


(- 1 ) 


fc-i 


1 / -k 
~ k \k- 1 
1 (2k - 2 
k\k— 1 J \ l — z 

Finally, we extract the coefficient of z n : 

1 (2k - 2 

1 (2k - 2 
k\k- 1 
1 (2k - 2 
k\k- 1 


z 


1-z 

k -1 


k -1 


‘ k ^ = i\k 


z 


n—k+l 


](1 


si -k 


1 -k 

n — k + l 
n — 1 
n — k + l 


- 1 ) 


n—k +1 


□ 


Combinatorial proof of Theorem \3.10 . It is well known that the num¬ 
ber of plane rooted trees (without multiple edges) with k vertices is 
given by the Catalan number Ck -i = \ Ck-f) ■ such tree can be 

transformed into a multi-edge tree of size n by distributing the n edges 
among the k — 1 edges of the non-multi-edge tree. This corresponds to 
a composition of n into k — 1 parts. There are (?~g) such compositions. 
Thus there are 1 ( 2 uZ\) (((Zy) pl ane rooted multi-edge tree of size n with 
k vertices. □ 
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The distribution of the number of vertices can now be derived from 


the explicit formula in Theorem 3.10 using Stirling’s formula. In order 
to determine the asymptotic behaviour of the moments, we use an 
approach via Hwang’s quasi power theorem which turns out to be more 
convenient. 

Theorem 3.11. Let V n be the number of vertices of a random rooted 
plane multi-edge tree of size n. Then 

4 9 /1 

E{V " ] = 5 n + 10 + ° (n 

V ™ = + + + 0 (1 


and 


P 


K 


4 

5 n 


n 


< v ) = 


J-c 


e -^/2 dt + 0 - L 


holds uniformly for v G M. Furthermore, the local limit theorem 


2y/2 


nil 


: expl — 


1 f k 


fn\ 2 

5 


h' n 


P(K = k) 

holds for k = + o(n 2//3 ). 

Proof. Let T(y,z ) be the bivariate generating function as in the first 
proof of Theorem 3.10 The functional equation (25) is equivalent to 

(26) zT(y, zf - (1 - z)T(y, z ) + y( 1 - z) = 0. 

Solving this quadratic equation yields 

(1 - z) - y/1 - z^l - (Ay + l)z_ 


(27) 


T(y,z) = 


2z 


note that the negative sign has to be chosen to obtain regularity at 
z — 0. 

The probability generating function of V n is then 


Pn(y) = 


[z n ]T(y,z ) 
[z n ]T(l, z) 


For y in a neighbourhood of 1, the dominant singularity is at z = 
1/(1 + Ay). As 


T{y,z) = 


1+4 y 


1+4 y 


a/1 - (Ay + 1 )z + 0(1 - (Ay + 1 )z) 


1+4 y 1+4 y 

= 2 y- y/y(l + Ay)y/l - (Ay + l)z + 0(1 - (Ay + 1 )z) 
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for z -A 1/(1 + Ay) except for one ray, singularity analysis [8j yields 
\z"}T(y, z) = + 1)"" _3/2 + 0((4j, + 1)“«- 5 ' 2 ) 


+ %) 

2a/ 


7rn J 


(% + l) r 


l + 0(- 
n 


For y — 1, this coincides with ({T|). 
Thus 


[z"]T(l,z) ^ 


4y + 1 


n+l/2 


1 + 0 ( - 
n 


The asymptotic formulae for mean and variance in Theorem |3.11| as well 
as the central limit theorem are an immediate consequence of Hwang’s 
quasi power theorem [12] in the version of [ITT, Theorem IX.8]. 

The local limit theorem follows immediately from the explicit formula 
in Theorem 3.10 applying Stirling’s formula to (23), we find that the 
total number of multi-edge trees with n edges and k = -jr + R vertices 
is equal to 


gn+3/2 


47TV^2' 


n z 


25 R 2 /l R R 3 \ 
exp i - —— + O ( - H-b — ) 


8 n 


n 


n 


n 


2 ) 


Combining this with the asymptotic formula 0 for the total number 
A n of multi-edge trees with n edges, we obtain the desired statement 
for R = o(n 2//3 ). □ 
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